Static charged fluid around a massive magnetic dipole 
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An analytical solution of Einstein-Maxwell equations with a static fluid as a source is presented. 
The spacetime is represented by the axially symmetric Weyl metric and the energy-momentum tensor 
describes a coupling of a fluid with an electromagnetic field. When appropriate limits are performed 
we recover the well-known solutions of Gutsunaev-Manko and Schwarzschild. Also, using Eckart's 
thermodynamics, we calculated the temperature, the mechanical pressure, the charge density and 
the energy density of the system. The analysis of thermodynamic quantities suggests that the 
solution can be used to represent a magnetized compact stellar object surrounded by a charged 
fluid. 
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I. INTRODUCTION 



The study of magnetic fields in astrophysical objects 
such as white dwarfs, neutron stars, pulsars and black 
holes, has grown sharply in recent years [H, [3, EJ Bl IH- 
In fact, several observations show that there are vari- 
ous scenarios where the magnetic fields and general rel- 
ativity can not be neglected . One of them is the pres- 
ence of strong magnetic fields in active galactic nuclei 
@) H H, • These nuclei are known to produce more radi- 
ation than the rest of the entire galaxy and directly affect 
its structure and evolution. Another scenario is the pro- 
duction of relativistic collimated jets in the inner regions 
of accretion discs, which can be explained considering 
magneto-centrifugal mechanisms [13, HH, EH O [3 [H| • 
Also, magnetic fields are important in understanding the 
interplay between magnetic and thermal processes for 
strongly magnetic neutron stars [1, [H, QjJ • At least 
10% of all neutron stars are born as magnetars, with 
magnetic fields above 10 14 G [H, OJ, Eij]- Analyti- 
cal models that describe these astrophysical objects are 
often associated with solutions of Einstein's equations 
[U S3, [H, [H, HE H||. In the context of relativis- 
tic hydrostatic models, the Tolman-Oppenhcimer-Volkov 
[27L l28j equations describe the internal structure of gen- 
eral relativistic static perfect fluid spheres, e.g. neutron 
stars. In the search for more realistic models for compact 
stellar systems, the energy-momentum tensor, the source 
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of Einstein's equations, is modified by introducing more 
complex terms that take into account additional physical 
properties as, for example, electromagnetic fields. In the 
case of relativistic magneto-hydrodynamics the reduction 
of this non-linear system leads often to simple models 
with limited applicability. For this reason, they are fre- 
quently superseded by numerical models [23, Ho, [H| • De- 
spite of this complexity, sometimes we can have useful 
analytical solutions, e.g., the Gutsunaev-Manko solution 
that describes the gravitational field of one static mas- 
sive magnetic dipole [32| . The aim of this work is to 
construct an analytical solution to the Einstein-Maxwell 
field equations coupled with a fluid in order to represent a 
static configuration that can be used to characterize the 
gravitational field of a magnetized astrophysical object 
surrounded by a charged fluid. 

The article is organized as follows. In Sec. |TT] 
we present the Einstein equations and the energy- 
momentum tensor to be considered. In Sec. Mil we ex- 
hibit the solution of the Einstein equations. In Sec. IIVI 
we study the thermodynamical properties of the system. 
Finally, in Sec. [Vj we summarize our results. 



II. COUPLING OF FIELDS 

The spacetime for our model is represented by the Weyl 
metric, 

ds 2 = f- 1 (e 2A [dp 2 + dz 2 ] + p 2 dip 2 ) - fdt 2 , (1) 

where / = f(p, z) and A = A(p, z). The coordinate range 
are the usual for axial symmetry. Our conventions are: 
G = c = 1, metric signature +2, partial and covariant 
derivatives with respect to the coordinate denoted 
by , /i and ; p, respectively. Greek indices run from 1 



2 



to 4, with (1,2,3,4) = (p,z,(p,t), and Latin indices run 
from 1 to 3. The aim of this work is to solve Einstein's 
equations with an energy-momentum tensor representing 
the coupling between a fluid and an electromagnetic field. 
Thus, in our model, the energy-momentum tensor is the 
sum of the electromagnetic energy- momentum tensor and 
a fluid energy-momentum tensor. The electromagnetic 
energy-momentum tensor considered is 



(EM) l 



-in 



F a pF^ .gapF^ijF 11 



(2) 



where F a p is the electromagnetic field tensor defined as 



F a g — A a 



1/3, a 



and is the four-potential = 



(U, V, W, $), where the functions U, V, W, $ depend only 
on the coordinates (p, z). The energy-momentum tensor 
of the fluid is 

T a p = £V a Vp + (p-<;9)h a p+T cl 0, (3) 

where v a represents the 4-velocity of the fluid, e is the 
fluid energy density, p is the fluid pressure, q is the bulk 
effective viscosity, 9 = v a a is the expansion, r a p is the 
stress tensor defined as r a p = —2i]a a p + q a vp + qpv a , q a 
is the heat flux, r] is the shear viscosity, h a p = v a vp +g a fj 
is the spatial projection tensor and a a p is the symmetric 
trace-free spatial shear tensor given by 



(4) 



The Einstein's equations for the system fluid plus elec- 
tromagnetic field are 



G a /3 = 87T 



(EM) 

T a p 



(F) 
T a /3 



(•5) 



(EM) (F) 

In the next section we solve Eq. with T a p and T a p 
given by @ and respectively. 



III. A SOLUTION TO THE EINSTEIN'S 
EQUATIONS 

To solve the Einstein equations we use a reference 
frame co-moving with the fluid. In this reference frame, 
the four-velocity of the fluid is v a = [0, 0, 0, V4 = V7]- 
Note also that 9 = v a a = in this frame, i.e. the ex- 
pansion of the fluid is null and there is no divergence or 
convergence of the fluid world lines. For this reason a 
co-moving observer does not see an effective spatial elec- 
tric current. Therefore, in our co-moving frame we have 
that the electric current is null, J 1 = 0. The condition 
J 1 = J 2 = 0, together with Maxwell equations, 



4ttJ^ = V Q i^ a , 



(6) 



leads to 



where hq is a constant, whereas the condition J 3 
be expressed by 



V • 



-vw 



= 0, 



can 



(8) 



where V = p-§^ + z-^;- Also, with the help of J 4 obtained 
from ((6|), we can write the charge density of the system 



-J^v^ as 



47rpe 



2A 



/ 



(9) 



Another simple equation can be obtained from the com- 
ponents {1,3} and {2,3} of Einstein's equations. Both 
components are equal to zero and can be written as 



(V p - U Z )W <Z = (V p - U z )W tP = 0. 



(10) 



One possible way to satisfy Eq. (|T0|) is to set the mag- 
netic potential W = constant. With this condition Eqs. 
|8]) and (flTjf are satisfied, but this leads to a physical 
model with tensions instead of pressures, and moreover 
these tensions are independent of the coordinate z. So, 
we discarded this solution. Another possibility is to set 
the magnetic potential W 7^ constant. With this last 
condition Eqs. (flU|) are satisfied when 



Vp - U- 



0. 



(11) 



Therefore, the constant ko in ([7|) is equal to zero. Now, if 
we add Einstein's equations {1,1} and {2,2} we find that 
the fluid pressure is equal to 



P = 



f 2 



47re 4A 



(v P -u z ). 



(12) 



from which we obtain, using Eq. (fTTj) , that it is equal 
to zero. From these considerations, the Einstein equa- 
tions reduce to two parts, one integrable system for the 
functions f,W,$> and A given by 



a =P.i£p Qi + L (w 2 - W 2 ' 



f 



K Z = P M^ + 2 Mp W, z - *, 
9 f 



A 



2/2 

A„.„ 



4 p 



= V • 



-VW 

p 



(13) 
(14) 

I 

(15) 
(16) 



pfe- 2A {V p - U lZ ) = k , 



(7) and another system of equations for the fluid density en- 
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ergy and flux radiation, which can be written as 



e = 



4irpe 2A 
qi = 0, 
Q2 = 0, 

v7 



<7:s 



47re 



2A 



(17) 

(18) 
(19) 

(20) 



where (3 is defined by 



P=l$- (7V/I - - VW ■ W - ^ V$ • V$. (21) 
2 \f J P J 

In the next subsection we find a solution for /, W, $ and 
A. 



A. Solution to the integrable system 



In order to solve the system of equations (I13H16|) . we 
first use Eqs. (fT3|) and (fl4|) to find the integrability con- 



dition A. p . z — A. 2 . p . This condition can be written as 



(22) 



Furthermore, if we substitute Eqs. (fT3"|) and (jT^J) into 
([TBI we obtain that 



^00* n ( P - 



(23) 



Assuming that the fluid energy density (e) and charge 
density of the system (<i c ) are different from zero, which 



implies from Eqs. © and (IT]) that /3 and V ■ ( | V$) 

are also different from zero, we find from Eqs. (|22|) and 
(HP that 



(24) 



which is satisfied when $ = <f> (/). Using this last condi- 
tion, we can rewrite equation (|22p as 



(2/*,/ - $) V ■ Q V<&) = i V • ^ V/ 

--VW • W - V • 1 ^$v$ 
p V/ 



(25) 



The explicit functional form of <&(/) is obtained in Ap- 
pendix A, replacing the function / by a power series of 
<f> in order to satisfy Eq. (|TB]). The result is $ = $ v7> 
where $0 is a constant. Now, substituting the functional 
form of $ into Eq. ([23)1 we obtain 



1 



VW- Vff = 0, 



(26) 



2 ' V/ / P 
which has the same form of the Gutsunaev-Manko's equa- 
tion for the massive magnetic dipole in vacuum if we de- 
fine W = yX — < J >2 -4, where A is the electro mag netic po- 
tential considered by Gutsunaev and Manko [32| • The so- 
lution of Eqs. (fT6|) and (f26|) written in prolate ellipsoidal 
coordinates x = (r + + r_)/2fc and y = (r + — r_)/2fc, 
where r± = \J p 2 + (z ± k) 2 and k is a real constant, is 
studied in the Appendix B. With this solution, the func- 
tions /, <&, W and A can be written in terms of only one 
parameter a as 



/ = 



4> 



W 



,2A 



1 / [x 2 — y 2 + a 2 (x 2 



l)] 2 + 4a 2 x 2 (l -y 2 )\ 2 



1 V [x 2 -y 2 +a 2 (x~ l) 2 ] 2 



4a 2 y 2 (x 2 



$0 



T / [.t 2 - y 2 + a 2 (:r 2 - 1)] + 4a 2 a; 2 (l - y 2 ) 
T I [.t 2 - y 2 + a 2 (a; - l) 2 ] 2 - 4a 2 y 2 (:r 2 - 1) 



4fca 3 x /T^$! (1 - V 2 ) [2 (a 2 + l) x 3 + (l - 3a 2 ) x 2 



y 2 + a 2 ] 



{a 2 + 1) {[x 2 ~y 2 + a 2 (x 2 - l)] 2 + Aa 2 x 2 (1 - y 2 )) 



(a; 2 - l) {[ x2 -y 2 + a2 ( x2 - !)] 2 + 4a 2 £ 2 (l - y 2 )) 



(a 2 + l) 8 



(x 2 — y 2 ) 9 



(i-*o) 



(27) 
(28) 
(29) 

(30) 



In the next subsection we relate the parameter a with B. Asymptotic solution 

the magnetic dipole of the system. 

Let us study the behavior of the metric and the elec- 
tromagnetic fields far from the source. In this case, it 
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is appropriate to write the functions /, $, W and A 
in spherical coordinates and expand them in power se- 
ries in r _1 . The spherical coordinates are related to 
the prolate spherical coordinates through the expressions 
x = {r—m)/k and y = cos#, where m is a real parameter. 
The function / far from the source takes the form 



/=1 



2fc(-l + 3a 2 ) 
(1 + a 2 )r 

2fc(-l + 3a 2 )(-fc + m + 3ka 2 + ma 2 ) 
(1 + a 2 ) 2 r 2 



+ 0(r~3). 
(31) 



Now, we imposed to the function / to have a 
Schwarzschild form far from the source, we find 



m(l + a 2 ) 
l-3a 2 



(32) 



Note that with this form of k, the third term in Eq. 
([3lTl is equal to zero, so finally we obtain that f = 1 — 
2m/r + (3(r~ 3 ). The function e 2A goes as e 2A = 1 + 
0(r~ 2 ). Therefore, far from the source we obtain the 
Schwarzschild metric. Hence m represents the mass of 
our system. 

Now, we analyze the asymptotic behavior of the elec- 
tromagnetic potential, — [0, 0, W, $]. For this pur- 
pose, we use the components of the four-potential (W, <&) 
written in spherical coordinates. Expanding Eqs. (|28|) 
and (129j) in series of r _1 we obtain that 



[r = B,/R^^ + Qni (33) 



$ = $o - 



(1 - 3a 2 ) 2 r 
^ + 0(r- 2 ). 



(34) 



Comparing Eq. (|33|) with the classical magnetic poten- 
tial, we note that the magnetic dipole moment is equal 
to 



(l-3a 2 ) 2 ' 



(35) 



which differs fro m the G utsunaev-Manko magnetic dipole 
by a factor of y/l — <&q. From Eq. ([55)) we can relate the 
parameter a with the magnetic dipole moment. Note 
that the parameter $o can take values between -1 and 
1. Therefore, from Eqs. (f33|) and (f34|) . we see that for 
values of |$o| near zero, the magnetic potential W and 
the electric potential $ far from the source attained their 
higher and lower values respectively. For values of |$q| 
near one, we have the opposite behavior. Note that the 
fluid total charge is m$o- 



IV. THERMODYNAMIC PROPERTIES 

A thermodynamic analysis of the system can be made 
performing a decomposition of the energy-momentum 



(EM) (F) 

tensor T a fj = T a g + T a p in its proper components as 
made by Eckart [33]. Note that for static system Eckart's 
thermodynamics does not present causality problems. 
Following Eckart, we can write the energy- momentum 
tensor in its proper components as 

T al3 = £,v a vP + Q a v + Q v a 



with £, Q a , \ a P are given by 

£ = v a vpT a P , 
Q a = -h«T^v 7 , 
X al3 = h*bP & T^. 



(36) 

(37) 
(38) 
(39) 



where £ is the energy density of the system, Q a is the 
energy flux of the system and is the stress tensor of 
the system. Calculating explicitly the components of the 
energy flux of the system, with the help of Eqs. (fT5)) and 
(|20|) . we found that Q 1 = Q 2 = Q 3 = which means that 
the system is in an equilibrium configuration. In this 
case the temperature of the system obeys the relation 



For a — 1, 2, 3 we have that 



0. 



(40) 



<9 P In 



d z ln 



TVf' 



= 0, 



(41) 

(42) 
(43) 

Using the above equations we state that the temperature 
of the system is equal to 

T 

J- DC: 



d v T = 0. 



T(fi,z) 



V7' 



(44) 



where is the temperature at infinity. 

The mechanical pressure of the system (P) is given in 
Eckart's thermodynamics by the expression 



487rpe 2A \f N J 



(45) 



where N = (2 - 3$g)/(2 - 2$g). The energy density of 
the system (|3"T1) can be written as 



16irpe 



(46) 



where M = (2 + 3$g)/(2 + 2$g). 

In Fig. [1] we present a graph of the temperature, the 
mechanical pressure, the charge density and the energy 
density of the system. Note that all variables decay 
rapidly to their asymptotic values, which suggest that we 
can treat our system as a compact stellar object. While 
varying the values of the parameters <&o, ot and m we 
found that we can obtain models which are less or more 
compact. For higher values of $0 we obtain the more 
compact objects. 
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Temperature Field Pressure 





V. CONCLUSIONS 

In this work we found an analytical axially symmet- 
ric static solution of Einstein's equations with a energy- 
momentum tensor which couples an electromagnetic field 
and a fluid field. Far from the source, our metric is consis- 
tent with the Schwarzschild solution. Moreover, if we let 
$o = we recover the Gutsunaev-Manko solution. The 
thermodynamic variables studied in our model suggest 
that we can treat our model as a compact stellar object, 
and these variables also may allow a direct comparison 
with observations. 

The solution to the Einstein equations presented in this 
work is a static metric. This looks like to be in contradic- 
tion to previous results that state that in the spacetime 
associated to a static magnetic dipole in the presence of 
electric charges appears a frame-dragging effect [H, [37| , 
i.e. we have a stationary metric. In our case, the static 
configuration is possible because the heat flux and the 
electromagnetic energy flux PD|) compensate each other 



in such a way that the total energy flux is equal to zero. 
Static metrics associated to the equilibrium configura- 
tions obtained from the superposition of counterrotating 
fluxes in the context of General Relativity are not new, 
e.g, in the Morgan & Morgan disks [38[ we have the same 
number of particles rotating in opposite directions (coun- 
terrotating hypothesis). 

We conclude by mentioning that our solution that rep- 
resents a charged fluid around a massive magnetic dipole 
may be useful to model a magnetar within a fluid. 
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APPENDIX A: FUNCTIONAL FORM OF $(/) 

To find the functional form of <I> first we write the 
function / as a power series of <J> 



(47) 



where a n are real coefficients. Replacing (j47|) in (fl"6|) . we 
obtain 



n=0 



■VVK 



+aiV 



= a V 



1 



(48) 



To simplify the notation we define the quantities B n 



B n = V • 



— vw 

p 



(49) 



with n = 0, 1,2,3, ...,oo, which satisfy the recurrence 
relation 



B n = 2$5„_ 1 - $ 2 £? n _ 2 . 



(50) 



With the help of Eq. ([50]) we can cast Eq. (J4SJ) into 
the form 



^ «nB n = Bi y^(q„ - 2a„ +2 )n$ n 

n=l 
oo 

+B 2 ^(n-2)a„$"- 2 = 0. (51) 



n=0 



ri=2 



To obtain one condition from Eq. ([ST]) to help us to 
find the functional form $(/), we demand that a n = 
2a n+ 2 f° r all n > 1. This allows us to write the above 
series only in terms of B2 , say 



00 



E2n + 1 2 
$ 
on 



n=0 
n=0 



= 0. 



(52) 



The above equation is satisfied if we set j3 2 =0, so 



B 2 = V ■ 



p 



0. 



(53) 



Finally, by direct comparison between Eq. (fB"3")) and 
Eq. (fTB|) we obtain that the functional form of <f> is 
*5 = &aV7, where $0 is a constant. Note that by 
D'Alembert criterion the series of Eq. (|47f converges 
when |$| < V2- 



APPENDIX B: SOLUTION OF EQS. © AND (f26|) 

Using the prolate ellipsoidal coordinates defined as 

p = ky/x 2 — lyl — y 2 , z — kxy, (54) 



where k is a real parameter, we can write Eqs. (|16[) and 
(l26l) in the form 



/ 



1 - y" 



■W tX 



0, 



(55) 



f k 2 1-y 2 

i-y\ 2(l-$ 2 ) Wf 



f 



-/,«- 



A; 2 x 1 — \ 



WA =0.(56) 



Now, we define the auxiliary functions 7? and G, so that 

H = i { lZy 2 ° } W and G = §±±/. With these definitions, 
Eq. (fST))) is trivially satisfied by a function I(x,y) if 



' — — — rG// , ■ 



x + 1 



G 



— [{l-y 2 )H\ y . 



(57) 
(58) 



(x + lf 

So, Eq. (fB"6"|) can be written in terms of G, H and I, say 



9 ? 



G t x 



2G 



(X + 1) I,xl,y 

G(x-1) 



(59) 



The next step is to assume that G, H and 7 are functions 
of the form 



N 
D' 



H=— and G=— , (60) 



where TV, I?, 5 and M are polynomials in x and y of the 
form, 



Jm.> Jr. 



N(x,y)= E W J 2/ m 



j=0 m=0 

3max Jmax 3 



j—0 rn—0 

D(x,y)=^ ^ 1 8 jm x'y m 

j=0 m=0 

M(x,y)= >T £ /ij^iT 

j=0 m=0 



(61) 
(62) 
(63) 
(64) 



with r/jm, crj m , 5j m , /ij OT being unknown coefficients to 
be determined. Substituting the relations (|6"0|) into Eqs. 
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(|57| and (|58|) we obtain the set of equations 

(x + l)(DN y - ND jV ) -(x- 1)(MS, X - SM, X ) = 0, 



(x + iy(DN x - ND, X )+M 2 



M 



0. 



(65) 



The system of Eqs (j6"5)) . together with Eq. can be 
solved in a tedious direct comparison of the coefficients 
of the polynomials involved. The solution depends on 
the number of coefficients for each polynomial, i.e. the 
value of jmax- Furthermore, we impose that the solution 
obtained be symmetric on the plane z = and also that 
has the Schwarzschild solution as a particular case. The 
first physical solution is obtained when j max = 4. In this 



case the coefficients can be written in terms of only one 
parameter, a. The solution is 

N = 8a 3 xy (x - 1) , 

D = [x 2 — y 2 + a 2 {x - l) 2 ] 2 - 4a 2 y 2 (x 2 - 1), 



M 



S = Aa 3 



y 2 + a 2 (x 2 -l)]+4a 2 x 2 (l-y 2 ) 



2x 6 



(1 - 3a 2 ) 
(a 2 + 1) 



y 2 + a 2 



(a 2 + 1) 



(66) 



Using the solution (J66J), the functions /, $ and W are 
known. The function A is calculated integrating Eqs. 
CE3 and (HU). 
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